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It is shown that a recent critique (arXiv:1210.1130 and arXiv:1211.4639 ) concerning the analytical 
solution of the Rabi model is unfounded. 
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It was demonstrated in (l| that the spectrum of the 
quantum Rabi Hamiltonian (oj = h = 1, g > 0), 



a^a + g<r x (a + ad) + Acr z 



(1) 



consists of two parts, the regular spectrum with energy 
eigenvalues £7„ + g 2 ^ IN and the exceptional spectrum 
(Juddian solutions) with Ef* c + g 2 G IN, which may oc- 
cur for special values of the model parameters g and A. 
The index ± denotes the parity of the regular eigenstate 
belonging to E^ . The exceptional states (if they are 
present) are all doubly degenerate with respect to parity. 
The parity invariance of the model (Hr commutes with 
e ma a cr z ) is instrumental to derive the function G±(x), 



G±(x) 



E 

71 = 



K n (x)(lT 



(2) 



where the K n (x) are known functions of g, A and x. 
G + (x) (G-(x)) determines the regular spectrum of ^ 
in the subspace with positive (negative) parity, because 
G± (x^ ) = for real x^ if and only if + g 2 = x^ . 

Maciejewski et al. argue in [2j that this result is in- 
valid because some of the real zeroes of G± (x) may not 
correspond to eigenvalues of (JT]). The authors of [2j do 
not dispute the fact that all points of the regular spec- 
trum correspond to zeroes of ([2]), but suspect that not 
all those zeroes are physical. I shall now prove that this 
is not the case. 

It is sufficient to confine the discussion to fixed (posi- 
tive) parity (negative parity is obtained by replacing A 
with — A in the subsequent formulas). The subspace 
with positive parity is isomorphic to B, the Bargmann 
space of analytic functions (see [l|) and the derivation of 
([2]) starts with the following system of coupled differen- 
tial equations for the wave function ip(z), which solves 
the Schrddinger equation H + ip(z) = Eip(z) in 7-L+, 



( z + 9)^Mz) + (gz ~ E)M*) + =0 (3a) 

(z - g)-^-fo(z) " & + E)Mz) + A&(*) = 0. (3b) 
dz 

Here we have used the notation yj(z) — 0i(z), ip(~z) = 
02 (z). This system corresponds to a linear homogeneous 




differential equation of the first order for the vector- 
valued function ^(z) = (0j (z), 2 (z)) T , 



(4) 



(j4|) has regular singular points at z = = ±17 and an irregular 
singular point (of rank 1) at infinity. Now it follows from 
the symmetry of this equation under the reflection z — > 
— z that the function $(z) = (02 (— z), 0i(— z)) T satisfies 
(j4]) as well. 01,2(2) have power series expansions around 

z = -9, 



h(z) = e-^J2K n (x)A 



{z + g) n 
x — n 



b 2 {z) = e-° z Y,K n (x){z + g y 



(5a) 
(5b) 



n=0 



where x denotes the spectral parameter x = E + g 2 . It 
follows that ^(z) is analytic in an open disk D\ with 
radius 2g centered at z = —g. Likewise, $(z) is analytic 
in a disk D2 with the same radius centered at z = g. 
All points in Do = D\ n D2 are ordinary points of (|4]) 
It means that if ^{zq) — $(zq) for any zq G Dq, 
^(z) and 3>(z) coincide for all z G Dq. But this entails 
that <f>\ (z) = ij)(z) is analytic in the whole complex plane 
because then 02 (~~ z ) is its analytic continuation beyond 
the radius of convergence of (|5a|) , comprising the second 
regular singular point of ((4]) at z = g. It follows that the 
conditions 



h(z ) 
h(z ) 



02 {-Zq) 
01 {-Zq) 



(6a) 
(6b) 



for any zq G -Do are necessary and sufficient for 4>(z) to 
be an element of the Bargmann space; the spectral pa- 
rameter x, determined by (|6a[l6bj) . corresponds therefore 
to an energy eigenvalue. Now (|6aj) is equivalent to (|6bj) if 
zq = G Dq, from which the expression for G±(x) given 
in Eq. @ follows immediately. This completes the proof 
sketched in In a comment [j| to the first version of 
this note, Maciejewski et al. still doubt the validity of the 
proof by invoking a standard theorem of complex analy- 
sis which says that a function holomorphic in a bounded, 



2 



connected region D of C vanishes identically if it vanishes 
at a denumerable infinity of points within D. This theo- 
rem has nothing to do with the present problem. Here we 
use the following elementary result 0] from the theory of 
linear differential equations: 

Theorem: Let the vector-valued function f(z) satisfy a 
linear homogeneous differential equation of the first order 
which has only ordinary points in the connected complex 
domain D. If f(z) vanishes at some point Zq G D, it 
vanishes everywhere in D. 

The condition f(zo) = &( z o) — ^( z o) — corresponds 
to (|6al I6bj) . and both are equivalent for zq = G Dq. 
If ([6a)) is satisfied at zq = 0, (|6b[) is satisfied as well 
and the vector / (0) vanishes. This is enough to con- 
clude that $(z) = ^{z) throughout Dq. It is not neces- 
sary that one of the components of f(z) vanishes at two 
distinct points (see below), but both components must 
vanish at one point. This is equivalent to the condition 
f(zo) — f'(zo) = if the scalar f(z) satisfies a second 
order differential equation, as the components of f{z) do. 

For z 7^ 0, Eqs. (|6a|) and (l6b|) are not equivalent and 
it becomes possible to have a solution to (|6a|) . while f|6b[) 
is not satisfied. This was discovered numerically in |2J for 
real Zq. Clearly, no unphysical solutions were obtained 
for zq = 0, but this is not a "lucky" accident as the 
authors of [2J believe, who checked the zeroes of G±(x) 
for x up to 30. 

The same argument applies to the generalized Rabi 
model with broken Z2-symmetry. Its Hamiltonian reads, 



oJa + ga x (a + a^) + ea x + Aa z 



(7) 



As was shown in |4j, the eigenvalue equation for H t 
is equivalent via integrable embedding to the follow- 
ing differential equation for the vector-valued function 

^{z) = (4>i{z),Mz)M*)Mz)) T , 



— = A(z)*(z), 



(8) 



with the coefficient matrix, 



-<>■■ 



A{z) 



\ 



z+g 




—A 

z-g 





E+e-gz 
z+g 
-A 
z— 9 






-A 

z+g 
E-t+gz 

z-g 





z+g \ 




E+e+gz 



(9) 



© has the same singularity structure as (Q| and re- 
gions Dq, Di, D2 can be defined as in the symmetric 
case. Due to the embedding, Eq. © has again a Z 2 - 
symmetry, which entails that with ^(z) also the function 
$(z) = (4>i(-z),(f>2(-z),<j)i(-z),(t)2(-z)) T satisfies ®. 
After expansion of *£(z) in powers of z around the reg- 
ular singular point —g, the condition ^(zq) = <&(zq) for 



zq G Dq leads to the following set of equations, 
AK~ 



-gz 



E- 

n=0 



E- 

n=0 ' 



Aif+ 



n=0 
00 

-(z + 9 )"=e^ 

n ^ — ^ 



ce-'-^K+izo+g)" 

n=0 
00 

Y j K-{zq+ 9 Y = 



-gzo 



n=0 



E 

00 

E 

n=0 



n=0 



AC 



a; — e — n 



AJf + 



*o) n , 
(10a) 

*o) n , 
(10b) 

(9-zqT, 
(10c) 

;(ff-*o) n , 
(lOd) 



with an unknown constant c. The are known func- 
tions of <7, A and x = E + g 2 . For zq = it is obvious 
that (gOll is equivalent to (fT0c|) and (|10b|) to (jlOdjl . We 
arc left with the two equations, 



E 

00 

E 

n=0 



K~ - 



x — e — n 
cA 



K- 



x ■ 



n 







.9" = 0. 



(11a) 
(lib) 



Eliminating c from Eqs. (11), we obtain the G-function 
for the generalized Rabi model [l[, 



with 



G e {x) = A 2 R + (x)R-(x) - R+(x)R-(x) 



R±(x)=J2Kt(x)g n 



R ± (x)=J2 



c — n ± e 
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(12) 



(13a) 



(13b) 



The reflection symmetry of the extended model allows 
to reduce the number of conditions as in the manifestly 
symmetric case. Therefore, the function W(x, g, A, e) de- 
rived in Q has exactly the same real zeroes as G e (x) and 
yields the same spectrum as seen in Fig. 5 of [2|. However, 
the proposed method is an interesting generalization of 
the approach introduced in [l[ which could be applicable 
to cases where embedding into a symmetric model is not 
possible. 

Regarding the numerical computation of the spectrum 
of the quantum Rabi model ([T]) , it may be advantageous 
to define a generalized G-function G±(x;z) by 



G±(x;z) = <fa{-z) - 4>i{z). 



(14) 



The vanishing of G±(x;zq) for zq G Dq corresponds to 
(pal) . Interestingly, this condition is sufficient to deter- 
mine the spectrum if 5(zo) 7^ 0. To see this, we note 



3 



that the conditions (6) correspond to a two-point bound- 
ary value problem for G±(x;z) in the complex plane, 
namely G±(x;z Q ) — G±(x; — z ) = 0. Because G±(x;z) 
satisfies a linear homogeneous differential equation of the 
second order, which is obtained from Eq. by elimi- 
nating ^2(2), this boundary value problem is incompat- 
ible [7| and has only the solution G±(x;z) = 0. Let 
G±(x; z) for j = 1,2 denote two linearly independent so- 
lutions of the differential equation satisfied by G±{x;z) 
(Eq. (11) in 0). Let us assume that G±(x;z ) = 
for z G Dq. Clearly, G±(x;z )* = as well. But 
because the coefficients of the power series of G±(x;z) 
in z are real (see Eq. (5)), we have G±(x;Zq) — 0, i.e. 
G±(x; z) vanishes at zq and Zq. This is again an incom- 
patible two-point boundary value problem if 3(zo) 7^ 0, 
because G±(x; z )G±(x; Zq) ^ G±(x; Zq)G\{x; z ) for 
almost all z G Dq. It docs not preclude isolated 
points zq G Dq which yield non-trivial solutions to 
G±(x;z ) — G±(x;zq) = 0. In [5| the following state- 
ment is made: "The key observation is that the func- 
tion in question satisfies a second order linear homoge- 
neous equation so that we only need to make it equal 
to zero at two distinct points." This is obviously in- 
correct in general; sin(z) satisfies the linear homoge- 
neous second order equation f"(z) — —f(z) and van- 
ishes at many distinct points without being identically 
zero: sin(z) solves the two-point boundary value prob- 
lem f(zi) = f(z 2 ) = e.g. for z\ — 0, z 2 = n. However, 
in the present case it is not required that G±(x;z) van- 
ishes at two points in -Do but that both components of 
f(z) — (G±(x; z), —G±(x; —z)) T vanish at a single point 
z G Dq. Because £ Do, the condition G±(x;0) = is 
necessary and sufficient for x — g 2 to be an eigenvalue of 
Hr. Maciejewski et al. write: " Numerical work seems 
to suggest that the condition at zero is somehow distin- 
guished, . . . " . The reason for this distinction is a simple 
mathematical fact explained above in great detail. 

Whereas isolated solutions of G±(x; z ) — G±(x; z^) = 
may exist for some zq G -Do, it cannot happen for 
zq G iK, because then zq — —Zq and G±(x;zq) = is 
equivalent to (|6a[ I6b[) . This is the interesting case from 
a numerical point of view, as it allows to overcome insta- 
bilities in the computation of G±(x;z) for large x. It is 
remarkable that G±(x; zq) has zeroes in x at the correct 
values even when zq G iM. is outside Do. The nonzero 
values of G±(x;zq) depend then on the order at which 
the defining series ([5a)) and (|5b| are truncated, but the 
position of x„ with G±(cc„ ; zo) = converges to the cor- 
rect value E^+g 2 for the following reason: The functions 
</>i(zo) and <p2(—zo) are holomorphic in C exactly at x„ , 
which entails a convergent series expansion in z at this 
point, even for z Dq. In Fig. [T]thc real and imaginary 
part of G+(x; 5?) is shown for x f» 71. zq is outside -Do but 
the joint zero of 5R(G+(a;; and 3(G+(a;; 5i)) allows to 



determine the energy eigenvalue — 70.00462935 to 
high precision, even though the series defining G+(x; 5i) 
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FIG. 1. Real and imaginary part of G+(x,5i) for g = 1 and 
A = 0.7 in the vicinity of x — 71. 

does not converge for x £ {a;^}- 

We conclude that the zq which may yield unphysical ze- 
roes of G±(x; Zq) are likely confined to the set IRn.Do\{0}. 
This set is of measure zero within -Do Q. In any case, 
its presence does not invalidate the results of [l[, which 
are based on the function G±(x) in Eq. ([2]) and not on 
the generalized function G±(x\z) for real z ^ 0, which 
is the subject of criticism in Q. The authors of @] use 
the Wronskian computed from <f>i(z) and §i{—z) together 
with their derivatives, which is equivalent with the condi- 
tions (6) and (10) for the Z2-invariant and the generalized 
Rabi model, respectively. The results of [2| neither cor- 
rect nor extend the findings of [l[ in these two cases, the 
only examples for which [2| presents explicit calculations. 
Nevertheless, Maciejewski et al. have correctly pointed 
out a gap in the derivation of G±(x) ]4^, giving me the 
opportunity to close it. 
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